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ABSTRACT 

Analysis  of  the  effect  of  multipath  propagation  on 
digital  communications  systems  was  conducted.  A  brief 
overview  of  the  root  causes  of  multipath  propagation  was 
included  in  this  discussion.  Probabilities  of  error  were 
then  derived  for  a  generalized  digital  communications  system 
experiencing  Rician  fading  due  to  multipath  propagation. 
Exact  results  were  obtained  for  equiprobable  M-ary  Frequency 
Shift  Keyed  and  M-ary  Phase  Shift  Keyed  modulation  schemes 
used  to  transmit  digital  data.  Furthermore,  the  probability 
of  error  was  obtained  for  a  generalized  digital 
communications  system  experiencing  single-bit  Intersymbol 
Interference  due  to  multipath  propagation.  Finally,  the 
performance  of  digital  communication  links  for  NASA's  Space 
Station  operating  in  the  presence  of  Intersymbol  Interference 
was  evaluated.  Results  obtained  tend  to  show  that  severe 
communication  system  performance  degradation  may  occur  on 
those  links  under  certain  transmitter/receiver  and  Space 
Station  geometries. 
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INTRODtJCTIOH 


This  thesis  was  intended  to  accomplish  four  purposes: 

(1)  To  introduce  the  reader  to  the  subject  of  multipath 
propagation  is  and  how  it  affects  the  performance  of 
digital  communications  systems, 

(2)  To  develop  as  general  an  expression  as  possible  for  the 
probability  of  error  of  a  digital  communications  system 
operating  in  a  Rician  Fading  environment  due  to 
multipath  propagation, 

(3)  To  develop  as  general  an  expression  as  possible  for  the 
probability  of  error  of  a  digital  communications  system 
that  experiences  Intersymbol  Interference  (ISI)  effects 
due  to  multipath  propagation,  and 

(4)  To  use  the  above  analysis  to  determine  whether  there 
would  be  significant  performance  degradation  due  to 
multipath  propagation  on  the  communication  links  for 
NASA's  proposed  Space  Station. 

In  a  certain  sense,  success  was  not  complete.  Pertinent 
points  follow. 

Chapter  2  does  indeed  contain  a  brief  introduction  to  the 
phenomenon  of  multipath  propagation.  As  a  matter  of 
practicality,  some  rather  important  concepts  were  kept  out  of 
the  description  in  order  to  enhance  brevity,  while  other  less 
important  items  were  included  in  order  to  make  the  analysis 
that  follows  in  later  chapters  more  understandable. 
Nevertheless,  Chapter  2  highlights  some  of  the  problems 
associated  with  communications  in  an  environment  where 
multipath  propagation  exists,  and  indicates  how  the  relevant 
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system  parameters  relate  to  the  physical  phenomena  associated 
with  multipath  propagation. 

Chapter  3  is  devoted  to  the  task  of  deriving  the 
probability  of  error  of  a  digital  communications  system 
experiencing  Rician  fading  due  to  multipath  propagation. 
Although  numerous  treatments  exist  of  generalized  digital 
communications  systems  experiencing  so-called  Rayleigh 
fading,  Rayleigh  fading  is  constrained  to  the  case  where  the 
reflecting  mechanism  is  dif f use--meaning  that  specular 
reflections  are  not  covered.  Since  multipath  propagation 
tends  to  occur  due  to  both  diffuse  and  specular  reflections, 
the  Rayleigh  Fading  model  is  in  essence  incomplete.  If  one 
is  to  consider  fading  caused  by  a  combination  of  both 
specular  and  diffuse  reflectors,  one  must  add  a  specular  term 
to  the  received  signal  model  resulting  in  a  Rician 
probability  density  function  (p.d.f.)  for  the  received  signal 
amplitude  statistics.  No  derivations  for  the  probability  of 
receiver  error  under  Rician  fading  have  been  found,  so  one 
was  attempted  here.  Rather  than  attempting  to  obtain  a  very 
general  result,  a  sub-case  of  significant  interest  was 
considered  in  detail.  Specifically,  only  equiprobable  signal 
transmission  was  considered,  and  a  complete  expression  was 
arrived  at  for  the  case  of  M-ary  Frequency  Shift  Keyed  (MFSK) 
modulated  signals.  Due  to  the  complexity  of  the  final 
expression  in  the  case  of  M-ary  Phase  Shift  Keyed  (MPSK) 


modulated  signals  which  prevented  us  from  showing  it  in  its 
full  form,  a  result  was  obtained  for  MPSK  as  well. 

Chapter  4  begins  with  a  presentation  of  the  configuration 
of  an  optimum  receiver  for  the  processing  of  one  of  M 
transmitted  signals.  The  probability  of  receiver  error  when 
ISI  effects  are  present  is  then  obtained.  In  order  to  keep 
the  derivation  manageable,  this  analysis  covers  only  the  case 
where  one  adjacent  bit  interferes  with  the  presently 
transmitted  bit.  Once  again,  the  final  result  covers  only 
the  case  where  the  M  possible  signals  have  an  equal 
probability  of  transmission.  Except  for  this  constraint,  the 
results  are  quite  general,  and  can  be  applied  to  a  larger 
class  of  modulation  schemes. 

Finally,  Chapter  5  addressed  the  problem  of  multipath 
propagation  as  it  applies  to  NASA's  Space  Station.  It 
includes  a  description  of  the  relevant  portions  of  the  Space 
Station  Communications  System  under  consideration.  Due  to 
the  fact  that  it  was  impractical  to  determine  appropriate 
fading  parameters  using  the  more  complete  Rician  model, 
multipath  analysis  for  Space  Station  was  carried  out  using 
only  the  ISI  model  developed  in  Chapter  4 .  Even  so,  and 
despite  the  fact  that  our  ISI  model  considered  only  the  case 
when  a  single  adjacent  bit  interferes  with  the  present  bit 
transmission,  it  was  demonstrated  that  there  may  indeed  be  a 
problem  with  poor  receiver  performance  under  certain 
conditions  when  the  Space  Station  is  operational.  This 
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Multipath  propigation  (or  just  multipath  for  short)  is  an 
insidious  perpetrator.  It  can  affect  the  performance  of  a 
communications  system  in  various  ways.  Trying  to  learn 
something  about  it  is  not  always  easy  because 
phenomenological  descriptions  require  knowledge  covering  a 
host  of  disciplines.  Physicists,  electromagnetic  engineers, 
and  communications  engineers  describe  it  differently  because 
all  three  tend  to  focus  on  different  aspects  of  this  physical 
phenomenon . 

Physicists,  not  surprisingly,  tend  to  focus  or.  the 
physical  structure  of  the  rut  ip  at  h  problem,  while 
electromagnetic  engineers  are  best  able  to  provide  relatively 
precise  descriptions  of  what  happens  to  the  phase  and 
amplitude  of  the  propagating  wave  for  a  set  transmitter  and 


an  EM  wave  from  some  point  "A"  to  some  point  "B"  and 
predicting  propagating  effects  is  often  too  complicated  to 
analyze.  Sometimes  one  may  find  that  the  only  way  to 
correctly  design  a  system  is  to  first  experiment  with  the 
channel  in  order  to  obtain  a  good  model  from  which  analysis 
and  design  predictions  can  be  obtained.  Nevertheless,  there 
are  fairly  general  theoretical  models  of  the  effects 
associated  with  multipath  that  communication  engineers  can 
use,  which  if  coupled  with  proper  understanding  of  the 
physics  of  the  situation,  may  result  in  reasonable  analyses 
and  performance  predictions.  The  intent  of  this  chapter  is 
to  provide  a  basic  introduction  to  the  phenomenon  of 
multipath,  with  particular  emphasis  on  the  salient  features 
which  are  pertinent  to  the  communication  problems  studied  in 
the  sequel. 

A.  CATEGORIES  AND  CLASSES 

As  described  by  the  physicist,  multipath  effects  are  a 
manifestation  of  either  reflection  or  refraction.  [1] 
Electromagnetic  engineers  also  include  diffraction, 
scattering,  focusing,  and  attenuation  as  modifiers  of 
multipath  propagation,  although  all  of  these  except  for 
diffraction  can  be  derived  from  the  basic  first  two.  [2,  p. 
2]  For  communication  engineers,  it  is  important  to  know  how 
those  mechanisms  affect  the  amplitude  and  phase  of  the 


received  signal,  as  both  can  influence  the  received  signal 
quality. 

Most  technical  sources  list  three  varieties  of 

electromagnetic  reflection:  specular,  resonant,  and  diffuse. 

From  the  communications  standpoint,  specular  and  resonant 

reflections  produce  the  same  results.  So,  for  our  purposes, 

it  is  sufficient  to  categorize  the  sources  of  multipath  as 

being  either  one  of  two  classes: 

•planar,  including  specular  and  resonant,  and 
•diffuse . 

Each  affects  the  communications  system  differently. 

1 .  Planar  R«f lacfcion 

This  is  the  term  used  for  cases  when  an  electromagnetic 
plane  wave  is  reflected  off  a  relatively  smooth  surface  by 
either  specular  or  resonant  (sometimes  called  "Bragg") 
scattering.  Planar  is  the  kind  of  reflection  most  sought 
after  in  applications  where  multipath  is  used  to  the 
advantage  of  the  designer,  such  as  when  ionospheric  "skips" 
are  used  in  order  to  extend  the  range  of  high  frequency 
transmissions,  since  it  usually  does  not  distort  the  signal 
significantly . 

Specular  reflection  is  what  occurs  when  the  surface  of 
the  reflector  is  very  smooth  in  relation  to  the  wavelength  of 
the  incident  wave.  Note  that  the  reflector  need  not  be 
smooth  at  all,  and  it  may  have  a  surface  that  on  the  whole 
appears  fairly  non-uniform.  However,  as  long  as  the  surface 
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appears  to  be  smooth  in  relation  to  the  wavelength  of  the 
incident  radiation,  (that  is,  surface  roughness,  p,  is  much 
less  than  the  incident  wavelength,  X)  specular  reflection  will 
result.  A  definition  of  roughness,  called  the  Rayleigh 
criterion,  is  as  follows 

4no  sin  (6) 

p  s  - —  (2-1) 

X 

where  O  is  defined  as  the  standard  deviation  of  the  surface 
irregularities  relative  to  the  mean  surface  height,  <|>  is  the 
angle  of  incidence  measured  from  the  grazing  angle,  and  X  is 
the  wavelength  of  the  carrier  [3]  .  Pictorially,  it  can  be 
viewed  as  shown  in  Figure  2-1. 

"Resonant"  reflection  occurs  when  a  surface  is  perhaps 
too  rough  to  reflect  specularly,  but  the  roughness  is 
periodic  with  a  spatial  separation  corresponding  to  some 
multiple  of  half  the  incident  wavelength.  When  this 
geometrical  scenario  is  set  up,  the  quanta  of  reflected 
energy  can  constructively  interfere  to  form  a  reflected  wave 
which  will  reach  the  receiver  [3] .  Resonant  reflection  can 
be  pictorially  viewed  as  shown  in  Figure  2-2. 

Again,  it  is  important  to  note  that  despite  the 
different  physical  phenomena  involved,  all  planar 
reflections,  both  specular  and  resonant,  affect  the 
transmitted  signal  in  basically  the  same  way. 
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generally  be  predicted  and  compensated  for,  unless  the 
structure  of  the  reflector  is  changing  with  time  in  an 
unpredictable  fashion.  In  such  a  case,  the  amplitude  change 
and  phase  shift  are  mathematically  treated  as  random 
processes,  resulting  in  a  single  ref lected-signal  path.  But, 
the  receiver  always  processes  a  single  direct-path  signal 
along  with  a  single  ref lected-path  signal  when  planar 
reflection  occurs. 

The  ref lected-path  signal  can  cause  a  phenomenon  known 
as  Intersymbol  Interference  (ISI) .  This  occurs  when  the  time 
delay  encountered  by  the  reflected  signal  is  so  long  that  it 
reaches  the  receiver  at  the  same  time  as  the  direct-path 
signal  for  subsequent  transmitted  symbols.  For  digital 
communications,  ISI  can  be  quite  destructive. 

Suppose,  for  example,  that  a  digital  communications  link 
set  up  with  bit  duration  T ,  releases  a  bit  from  the 

transmitter  at  time  t,  and  via  the  direct  path,  reaches  the 
receiver  at  time  t+8.  The  ref lected-path  signal  will  travel 

a  longer  route  to  get  to  the  receiver,  so  the  bit  traveling 
this  path  will  reach  the  receiver  at  a  later  time  t+8+ET.  If 

bit  "k"  is  transmitted  at  time  t,  the  subsequent  bit  "k  +  1" 
can  reach  the  receiver  via  the  direct  path  before  or  during 
the  time  that  the  reflected  bit  "k"  does.  ISI  will  then 
result.  This  situation  can  be  graphically  viewed  as  shown  in 
Figure  2-3 . 
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FIGURE  2-3.  Intarsymbol  Interference  (ISI) 

If  there  is  more  than  one  reflected  signal  path  causing 
several  bits  to  be  received  simultaneously,  the  situation  can 
become  almost  untenable.  This  will  be  a  topic  of  further 
discussion  in  Chapter  4  of  this  thesis.  Fortunately,  the 
study  of  ISI  has  received  significant  attention  in  the  past 
so  that  relatively  effective  techniques  to  combat  ISI  are 
available . 

It  is  important  to  note  that  if  there  is  a  single  planar 
reflector  present,  and  if  enough  information  about  the 
reflector's  characteristics  is  at  hand,  an  optimum  or  near 
optimum  signal  processing  receiver  structure  is  usually 
obtainable . 

2  •  Diifuaa _ Rmflmetion 

Multipath  propagation  due  to  diffuse  reflection  usually 
occurs  when  the  reflector  is  not  a  discrete  object.  With 
diffuse  reflection,  the  amplitude  and  phase  of  the  signal  can 
be  altered,  but  there  are  further  complications.  The 


received  signal  may  in  fact  be  a  continuum  of  reflections, 
with  its  incoming  phase  smeared  to  a  point  where  it  becomes 
undecipherable.  The  amplitudes  of  any  one  reflected 
component  may  be  indistinguishable  from  the  others.  So- 
called  Rayleigh  fading  is  a  common  manifestation  of  diffuse 
reflection.  The  ionoshpere,  for  example,  may  act  as  a 
diffuse  reflector  to  certain  signal  frequency  bands.  Some 
troposhperic  effects  such  as  precipitation  and  ducting  can 
behave  as  diffuse  reflectors.  Thus,  although  diffuse 
reflection  is  not  a  simple  phenomenon,  it  is  common  and  must 
be  dealt  with.  Fortunately,  this  problem  received  early 
recognition  by  analysts  and  system  designers  so  that  the 
available  pertinent  analytical  tools  have  reached  a  fairly 
mature  stage. 

Since  diffuse  reflectors  seem  to  sometimes  act  as  a 
continuum  of  individual  specular  reflectors,  computer  models 
can  be  developed  in  which  the  diffuse  reflector  is  simulated 
as  a  sum  of  a  large  number  of  specular  reflectors. 

As  with  specular  reflection,  the  characteristics  of 
diffuse  reflectors  can  change  over  time  and  appear  to  be 
random,  thus  necessitating  the  use  of  random  processes  in  the 
mathematical  model. 

B .  FREQUENCY  EFFECTS 

The  physics  of  multipath  are  fairly  straightforward.  If 
the  reflection  occurs  off  a  fairly  flat  plate,  the  effect 


will  appear  to  be  specular.  If  the  reflector  is  less 
discrete,  diffuse  multipath  effects  can  be  expected. 

But  what  appears  flat  and  smooth  to  one  packet  of  EM 
waves,  may  appear  differently  to  another.  For  instance,  if 
long-wavelength  high  frequency  radiation  were  bounced  off  a 
large,  flat  structure  that  appeared  smooth,  the  radiaticr. 
would  be  specularly  reflected.  However  the  surface  roughness 
of  the  structure  may  be  such  that  it  might  appear  rough  tc 
shorter-wavelength  EM  waves  and  cause  diffuse  reflection. 

In  reality,  no  reflector  is  either  purely  planar  or 
diffuse.  But  in  many  instances,  the  problem  can  be 
simplified  to  one  in  which,  for  a  given  set  of  frequencies, 
a  reflector  can  be  modeled  as  either  planar  or  diffuse 
because  one  type  of  reflection  is  clearly  dominant. 

C.  MULTIPATH  IITHCTH 

There  are  two  basic  ways  to  describe  the  effects  of 
multipath  from  a  communications  standpoint.  They  are:  (a) 
Fading,  and  (b)  Intersymbol  Interference  (ISI). 

When  a  transmitted  signal  arrives  at  the  receiver,  the 
transmitted  bit  is  adjoined  by  scaled  and  delayed  versions  of 
the  transmitted  bit  due  to  multipath.  These  scaled  ana 
delayed  versions  will  interfere  with  the  normal  performance 
of  the  receiver.  It  is  possible  to  empirically  measure  the 
degradation  in  receiver  performance  under  given  multipath 
conditions  in  terms  of  signal  loss,  usually  expressed  :r. 


units  of  decibels.  When  the  receiver's  measured  signal  loss 
is  expressed  in  this  fashion,  the  effect  is  called  multipath 
"fading." 

Another  common  method  of  determining  the  effect  of 
multipath  on  a  receiver's  performance  is  to  analytically 
express  the  probability  of  error  for  the  receiver  under 
circumstances  where  a  previously  transmitted  bit  is  received 
at  the  same  time  as  the  presently  transmitted  bit.  That  is, 
adjacent  bits  can  be  delayed  by  the  multipath  to  such  an 
extent  that  they  are  received  at  the  same  time  as  bits  that 
were  transmitted  later.  This  phenomenon  is  called 
Intersymbol  Interference,  or  ISI.  This  method  is  purely 
analytical,  and  therefore  can  be  insoluble  for  certain 
conditions  of  multipath  whereby  a  large  number  of  bits  are 
delayed  in  such  a  fashion  that  many  such  bits  are  received 
simultaneously . 

It  is  possible  to  rule  out  any  form  of  multipath 
degradation  if  the  following  two  criteria  are  satisfied: 

(1)  There  are  no  RF  reflectors,  either  specular  or  diffuse, 
which  generate  multipath  between  the  transmitting  and 
receiving  platforms,  and 

(2)  The  receiving  platform  is  small  enough,  simple  enough, 
or  "clean"  enough  so  that  the  transmitted  signal  is  not 
"bounced  around"  within  the  platform's  structural 
boundaries  before  the  signal  reaches  the  receiving 
antenna.  It  is  important  to  note  that  this  criterion 
is  heavily  dependent  on  the  bit  rate,  because  the  bit 
rate  is  what  determines  whether  the  receiving 
platform's  structure  is  of  a  size  that  will  result  in 
this  multiple  bounce  effect  described  above. 
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D.  SIGNAL  MODKL 


The  general  mathematical  form  of  an  M-ary  digital 
communication  signal  is 

s .  T(t)  =  J2  f.(t)  C0S{(Dct  +  f(t)}  (2-2) 


where  the  subscript  i  represents  any  one  of  M  possible 

transmitted  signals,  the  subscript  T  represents  the  duration 
of  the  signal,  f^(t)  determines  the  amplitude  shape  of  the 

i'th  signal,  C0c  is  the  carrier  frequency,  and  <j>i(t)  is  the 
deterministic  phase  of  the  signal.  This  mathematical  form  is 
equivalent  to 


s 


(t) 


Re 


rtt) 


j[0)ct+$.  (t)  ] 

e 


(2-3) 


Since  this  form  is  mathematically  easier  to  manipulate,  it  is 
sometimes  preferred  over  the  cosine  form  of  Equation  2-2. 

To  illustrate  how  these  expressions  are  used,  consider 
Binary  Phase  Shift  Keyed  (BPSK)  modulation.  In  this 
modulation  scheme,  there  are  two  possible  signals  to  be 
transmitted.  These  signals  would  be 

s,  T(t)  =  Jl  A  cos(G)„t  +  Jt)  (2-4) 

s;  _(t)  =  Jl  A  cos(U>„t  +  n)  (2-5) 

Thus,  for  this  modulation  scheme,  f^t)  would  in  both  cases 
be  the  signal  amplitude.  A,  and  signal  ”1"  is  differentiated 
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from  signal  "2"  by  the  opposite  (antipodal,  in  this  case) 
phases . 


Frequency  Shift  Keyed  (FSK)  modulation  is  similar,  but 
instead  of  substituting  a  fixed  phase  term  in  the  place  of 
4>i(t),  a  function  corresponding  to  the  integral  of  the 
instantaneous  phase  must  be  used. 

In  the  case  of  an  Amplitude  Shift  Keyed  (ASK)  modulation, 
the  <J>i(t)  terms  would  be  constant  or  eliminated,  and  the  f^(t) 

term  would  be  used  to  define  the  different  possible  signal 
amplitudes . 


E .  MULTIPATH  PARAMETERS 
1  .  Fading  Modal 

Chapter  3  presents  the  analysis  pertaining  to  the 
fading  effects  of  multipath  propagation.  Utilization  of  this 
analysis  requires  that  the  following  steps  be  taken: 

(1)  Experimentation  or  computer  modeling  of  the  reflecting 
structure  must  be  performed  in  order  to  provide  the 
necessary  multipath  parameters.  In  the  case  of  the 
Rician  model  presented  in  Chapter  3,  the  parameters  are: 

a  s  amplitude  of  the  planar  (specular)  component 
5  =  phase  of  the  planar  (specular)  component 
v  =  amplitude  of  the  diffuse  component 
0  =  phase  of  the  diffuse  component 
The  amplitude  and  phase  of  the  specular  component  can 
normally  be  considered  to  be  deterministic,  whereas  the 
parameters  pertaining  to  the  diffuse  component  are 
usually  considered  to  be  random  processes.  In  this 
analysis,  those  random  processes  are  assumed  to  be 
Gaussian,  since  this  is  the  case  that  most  often  occurs. 

(2)  Once  these  parameters  have  been  determined,  the 
modulation  scheme  to  be  used  must  be  substituted  for  the 
general  signal  model  used  in  the  analysis  to  follow. 


(3)  At  this  point  the  receiver  performance  may  be  obtained 
directly  from  the  results  and  the  probability  of  error 
curves  may  then  be  constructed. 

We  have  provided  the  probability  of  error  for  the  M- 
ary  FSK  (MFSK)  case,  as  well  as  the  steps  to  be  taken  to 
develop  the  M-ary  PSK  (MPSK)  results. 

2.  131  Modal 

Chapter  4  analyzes  the  effects  of  multipath  that 
result  in  ISI.  Although  it  may  be  possible  to  develop  a 
general  solution  for  performance  affects  due  to  multipath- 
induced  ISI,  the  problem  rapidly  becomes  unwieldy.  In  order 
to  keep  the  problem  manageable,  an  a  priori  assumption  was 
made  that  the  modulation  scheme  to  be  used  would  be 
orthogonal.  Therefore,  this  analysis  is  only  valid  if  an 
orthogonal  modulation  scheme  is  used. 

The  required  parameters  necessary  to  evaluate  the 
effects  of  ISI  are: 

Pi  s  amplitudes  of  the  multipath  signals 
"H i  =  travel  delays  of  the  multipath  signals 

These  parameters  can  be  obtained  analytically  through  simple 
geometric  approaches. 

The  ISI  analysis  will  be  illustrated  with  an  example 
in  Chapter  5. 

F.  LIMITATIONS 

Analysis  involving  ISI  effects  has  been  performed 
specifically  for  communication  links  to  be  used  in  NASA' s 

d? 


proposed  Space  Station.  The  results  are  presented  in 
Chapter  5 . 

A  specific  example  involving  fading  effects  is  not 
included  here,  since  to  utilize  the  performance  results 
obtained  in  Chapter  3,  an  empirical  model  of  the  physical 
configuration  must  be  available,  either  through 
experimentation  or  computer  modeling  of  the  reflecting 
structures.  This  can  be  (and  typically  is)  an  expensive  and 
time-consuming  process.  Nevertheless,  such  efforts  must  be 
carried  out  in  order  to  quantify  the  multipath  fading 
parameters.  This  fact  limits  somewhat  the  utility  of  the 
fading  performance  analysis  (excluding  its  purely  theoretical 
value)  because  the  lack  of  procurement  of  realistic  values 
for  the  fading  parameters  prevents  the  application  to  a  Space 
Station-related  problem.  An  analysis  related  to  the  Space 
Station  geometry  using  arbitrary  values  for  the  necessary 
parameters  would  be  of  little  value. 

Finally,  despite  the  fact  that  a  theoretical  analysis  of 
the  physical  structure  of  the  reflector  can  provide  useful 
parameters  for  use  in  the  analysis  of  ISI  performance 
effects,  since  the  modulation  scheme  considered  for  the 
communication  link  is  orthogonal,  the  analysis  and  results  of 
Chapter  4  do  apply. 
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III. 


A.  INTRODUCTION 

In  this  chapter,  we  carry  out  analysis  which  models  the 
effects  of  multipath  in  the  form  of  a  system  that  induces 
fading  on  signals  transmitted  over  such  a  communications 
channel.  At  first,  we  will  assume  no  particular  modulation 


scheme  in  order  to  ensure  that 


;he  analysis  w:l.  ce 


applicable  to  as  many  common i ^ at  ions  systems  as  possible. 
Thus,  a  mathematical  model  for  a  generalized  digital 


communications  signal  will  be  used  here. 


-  cl  0  r  w$  w  *  *  . 


conclude  the  analysis  by  specifying  the  resultant  probability 
of  error  { P { e  > )  for  a  system  employing  M-ary  Frequency  Shift 
Keyed  (MFSK)  modulation  as  well  as  the  mathematical 
expressions  that  lead  to  P{e}  for  M-ary  Phase  Shift  Keyed 
(MPSK)  modulation.  We  will  also  assume  that  no  multipath 
countermeasures  (such  as  diversity  techniques)  are 
implemented . 


B.  PRELIMINARY  ANALYSIS 


Assuming  no  particular  digital  modulation 


rheme , 


ae.neral  mathematical  form  of  the  transmitted  sicn.al  is 


s  .  ( t ) 


-yr 


(t)  e 


A* 


(3-1) 
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s  in  6 
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V  cos  0 


FIGURE  3-1.  Rician  Fading  Signal  Geometry 


Using  the  larger  triangle  in  Figure  3-1,  we  have 
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v1  =  a  +  v  +  2av  cos  (0-5) 


In  most  analyses,  it  can  be  assumed  that  a  and  8  are 
deterministic  parameters  while  v  and  0  are  random  variables 


with  joint  probability  density  function  (p.d.f.) 


fve(v,8)  =  *^7  exp  -  ~ 

V,e  2d  202 


(3-5) 


where  0<v^°°,  and  0^0^271,  so  that  v  and  0  can  be  seen  to  be 
independent  Rayleigh  and  uniformly  distributed  random 
variables  respectively. 

It  is  shown  by  Turin  [4}  that  due  to  the  given  joint 
p.d.f.  of  v  and  0,  the  joint  p.d.f  of  v'  and0'  is  given  by 


f  (v\  0')  =  - 

v'-0,  2no2 


exp  I- 


2  2 

v'  +  a  -  2av'  cos(0'-8) 


(3-6) 


where  0  <,  v'  <i  <»,  and  0  S@'-  8  <  2 K,  so  that  by  direct 
integration 


,<v. 


2  2 

v'  v*  +  a  -  2av  cos(0'-8) 

- -  exp - - -  d0' 

27KJ  20 


=  --  exp 


.  2  2 

v'  +  a 


a2  I  2a2  I  °  U 


(3-7) 


where,  0^v'<,oof  resulting  in  a  Rician  density  for  the  random 


variable  V' 


Whereas©  is  a  uniform  random  variable,  0'  is  not  uniform, 


as  can  be  seen  from 


£  f  v'  V 2+  a2-  2av’  cos  (0'-  5) 

f«,(9  )  -  - r  exp - - - -  dv’ 

e  2  2m2  2o! 


(3-8) 


Several  steps  must  be  carried  out  in  order  to  evaluate 


the  integral,  resulting  in 


x  -CL/2  O2 

V'9''  -  *  e 


ojln 


(Ql  E,  o2sin2{0,-8)  acos(0'-8) 

(0-0)  exp - - -  erf 4  - 

2o  o 


(3-9) 


where  0  <!  0 '  -  8  <;  271.  The  actual  shape  of  the  p.d.f 


dependent  on  the  CL/O  ratio.  Note  that  for  a  purely  diffuse 


multipath  contribution,  resulting  in  a/O  =  0, 


f  <0’>  =  T“ 
»•  271 


(3-10) 


i  •  ii  f  r  r 


which  is  a  uniform  p.d.f.,  where  0  £0'-  8  <  27t.  For  a  purely 
specular  multipath  component,  which  results  in  Cx/C— » 


f 

«• 


(0*)  - 


0'- 8  *  0 
0  5  =  0 


(3-11) 


which  corresponds  to  a  Dirac  delta  function.  A  complete 
analysis  would  require  checking  that  the  relationship 


(3-12) 


is  indeed  satisfied,  however  some  integration  difficulties 
arise.  Nevertheless,  this  clearly  demonstrates  that  9'  does 
not  have  a  simple  uniform  p.d.f.  (unless  a/a  =  0). 

Later,  the  case  of  8  being  a  random  variable  will  be 
investigated.  This  analysis  may  be  useful  if  the  specular 
reflector's  position  is  unknown  or  if  it  is  moving.  It  will 
thus  be  assumed  that  8  has  p.d.f. 


m  cos  (8) 

e 

f  (8)  =  - 

a  2711  (m) 


(3-13) 


where  — 7C<S<7C,  and  O^m^®®,  which  due  to  the  variable  parameter 

m,  represents  a  family  of  p.d.f. 's  ranging  from  the  uniform 
p.d.f.  for  m=0  to  a  Dirac  delta  function  for  m— »®®.  It  will 

be  necessary  to  determine  the  overall  effect  on  the  p.d.f.'s 


involving  v*  and  0 '  when  5  is  random.  From  previous  work, 
however,  even  if  5  is  random,  v'  will  have  an  (unchanged) 
Rician  p.d.f.  as  given  by  Equation  3-7. 

Returning  now  to  the  received  signal  s^ft),  we  have 


jco  t  j<j>.  (t) 

18  18  c  1 

<ae  +  ve  )  f .  (t)  e  e 

/|  j  (0)  t+$.  (t)  +8)  j  (CD  t+4>.  (t)  +0  ) 

2  Rel  af  .  (t )  e  +y2Revf.(t)e 

3.(t)  =  VT  a  f^t)  cosToj^c  +  (t )  +  sj 

♦VT  v  f.(t)  cos^(0ct  +  ^(t)  +  e]  (3-14) 


The  term 


yr  a  f.  (t)  cos^coct  +  «)>i  (t)  +  sj 


(3-15) 


is  the  so-called  specular  (multipath)  component,  but  as 
discussed  in  Chapter  2  of  this  thesis,  a  better  name  for  it 
would  be  the  "planar"  reflection  term,  since  it  can  be  caused 
by  either  specular  or  resonant  reflection.  In  all  our 
analyses,  a  will  be  treated  as  a  deterministic  quantity.  The 


term 


yr  v  f.(t)  cos|co.t  +  ()>.  (t)  +  sj 


(3-16) 


is  the  so-called  diffuse  (multipath)  component,  where  v  and  0 
are  random  variables  with  the  previously  specified  p.d.f. 's. 


Observe  that 
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s,(t)  =  [a  +  v  cos  (8-0)  ]  J 2  fi 

(t)  COS  [coct  +  <(>,  (t)  +8]  M 

+  v  sin  (8-0)  ■JT'  f.  (t)  sin[(fl 

ct  +  <t>;  (t)  +  5]  (3-17)  « 

For  convenience  we  define 

1 

a  =  a  +  v  cos  (8-0) 

(3-18,  1 

a2  =  v  sin (5-0) 

(3-19)  I 

Clearly,  ax  and  a2  are  random 

variables,  with  p.d.f.  9 

obtainable  from 

v  =  1  (a  -  a)  +  a. 

(3-20)  a 

8-0  =  tan  1 - 

(3-21) 

la  -  a 

-*1 

ana 

H 

9ai 

9 

-  =  cos  (8-0) 

dv 

(3-22) 

l 

dai 

Em 

~ —  =  v  sin (5-0) 

90 

(3-23)  " 

-  =  sin(5-0) 

(3-24) 

9v 

>« 

i 

9a2 

• , 

-  =  -v  sin (8-0) 

90 

(3-25)  £ 

••• 

The  resulting  Jacobian  of  the  trans: 

J 

formation  J(v,0)  is  given 

mm 

by 

§ 

'  i 
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% 

y 

> 

fi 

71 

V 

v‘  Q.  V  w . t-". A.1  *  1  >jiib  >»■!. 

J(v,0)  = 


cos  (8-0) 

v  sin(5-0) 

sin{8-0) 

-v  cos  (8-0) 

{  [<ara|!+  *,]  - 


5- tan 


2  2 

(aj-a)  +  a2 


2  2 

(a^a)  +  a2 


^\j  a  \  f  (a,~tx)  +  a  1  .  8-tan  k 


/ (a,-  a) 


U  (at-a)  +  a2J  ,  8-tan  ‘^a2/(a,-a)  J 


GJ2K 


2  2 

1  (ax-  a)  +  a2 

exp  17 


l  -(a1-a)Z/2a2  -a2/202 

— / - -  e  — ; - -  e  (3-29) 

<V2jt  oy27i 


This  demonstrates  that  a!  and  a2  are  independent  Gaussian 
random  variables  of  equal  variance,  G2 ,  and  means,  a  and  0 
respectively . 

Thus,  if  we  define 


yn(t)  s  /7  f  (t)  coskt  +  f(t)  +  sj 


(3-30) 


yi2(t) 


J2  f.(t)  si  r|a)ct  +  <(>.  (t)  +  s] 


we  have 


s,  (t) 


-t 


\  y*(t) 


(3-31) 


(3-32) 


where  i  =  and  Tb  S  t  ^  Tf 

We  now  define  some  quantities  of  interest,  namely 

J  y*,  (t)  dt  =  J  2T(t)  J-J  +  ~  cos  (20). t  +  20.  (t)  +  28)  J  dt 


f  (t)  dt  h  E,. 


(3-33) 


since  for  any  fi(t.)  that  varies  much  more  slowly  than  the 
f  requency  C0c 


J  f  ( t )  cos|2(0)ct  +  0.  (t)  +  8)  j  dt  *  0 


(3-34) 


Also, 


yi2(t)  dt 


'  J  [l  -  7 


ccs(2(0.t  +  20.  (t)  +  25)1  dt 


1 


(t)  dt  s  t,, 


(3-35) 
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Finally, 


J  yu  (t)  yi2(t)  dt 

T 

*b 


‘  f 

|2f2(t)  cos  (0)  t+0.  (t)  +  8)  sin  (CO  t+<j>  (t)  +  8)  dt 

J  1  Cl  C  '  1 

T 
*  b 


l 

=  J  f.  (t)  s  inj  2  (COct  +  4>.(t)  +  5)J  dt  = 


(3- 


‘b 


We  also  define 


P 


i  j 


i 

1 


ft(t)  f.(t) 


dt 


=  1 


(3 


Later,  we  will  encounter 


Tf 

1  yn(t)  dt 


=  J  2f.(t)fm(t)  cos  [<o,t+  0  (t )  +  8]  cos  [co  t+ 0_  (t )  +  5] 


T. 


‘J 


=  }f.(t)f  (t)  cos  [0.  (t)  -  0  (t))dt 


+ 


Jf,(t)ffl(t) 


cos[2o)t  +  0  ( t )  +  0  (t)  +25] 

c  i  fn 


dt 


(3 


36) 


37) 


dt 


-33) 
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The  second  integral  in  this  expression  is  clearly 
insignificant,  whereas,  the  first  integral  involves 


jf  (t> 


~j0  (t) 


Re  I  If,  (t)  e 


f  (t)  e 

JT\ 


dt  s  5,. 


(3-39) 


which  may  not  be  zero  even  for  i*m.  Furthermore, 


jyi2(t)  yr2(t)  dt 


(t)f  (t )  cos  [0.  (t )  -  0_  (t )  ]  dt 


-  [f  (t)f  (t)  cos  (20)  t  +  0  (t )  +  0  (t)+28)  dt  =  £  (3-40) 

I  1  W  C  1  m  1.T 


Also, 


y*2{t)  dt 


=  f;(t)f  (t)sin[0  ( t )  -  0  <  t )  ]  dt 


+  Jf ;  (t)  (t )  sin  [2(i)_t  +  0  (t )  +  <|Mt )  +  25]  dt  « 


(3-41) 
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where 


Jf^UfJUsin^Jt) -<|>.<t)]  dt  s  p.n 


(3-42) 


C.  THE  DETECTION  PROBLEM 

In  the  M-ary  hypothesis  testing  problem,  the  received 


signal  r(t)  is  modeled  by 


r (t )  =  s.(t)  +  n(t) 


(3-43) 


where  i  =  1,2,...,M,  Tb£t^Tf,  and  n(t)  is  the  Additive  White 

Gaussian  Noise  (AWGN)  of  Power  Spectral  Density  (PSD)  level 
N0/2  Watts/Hz.  If  s^(t)  is  a  completely  known  signal  with 
probability  P^  of  being  transmitted,  then  the  optimum 

receiver  must  compute 


In (P. : 


i. 

i 


exp 


-f 

N„  J 


r (t)  st(t)  dt 


-K 


(t)  dt 


(3-44) 


for  i=l,2,...,M,  and  will  decide  that  sm(t)  was  transmitted 
if 


1  >  1.  V  i  t  m 

m  l 


(3-45) 


Since  due  to  multipath,  the  received  signal  model  is 


2 

si(t>  =  X  a*  yik(t> 

k-l 


(3-46) 
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for  i 


=  1,2,... ,M,  and  Tb<t£Tf,  then  conditioned  on  the  random 
variables  ak,  k=l,2, 


ln(Pi>  [2  f  if 

Va2  =  e  exp  —  2A  Rik-  —  2/fi  a< 

4  0  0  k=l 


(3-47) 


where 


Ru  s  r(t>  yi*(t)  dt 


(3-48) 


and  the  simplifying  condition  (see  Equations  3-38  and  3-40) 


[  y..  (t)  y.  (t)  dt  = 

I  J  ik  m 


0  n  *  k 


E  n  =  k 


(3-49) 


has  been  used  in  Equation  3-48. 

Using  the  previously  derived  joint  p.d.f.  for  ax  and  a2, 


00 

xi  =  JJ  fA1,A2(ai'a2)  dai  d 

-00 

ln<P‘’  f  J  ,  2 

.  e  Jexpl-  —  a,  ♦  -  R:ia: 

[J  sn  i  2  ‘A  da.- 

J  T  N  2  N  -2  2  7  2  / 

—  L  0  0  20  J  V 


(a.  -  a)  I  da. 


(3-50) 


We  work  through  the  more  general  (first)  integral  of 
Equation  3—50.  The  exponent,  only,  is  of  the  form 
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Thus,  the  first  integral  is  of  the  form 


-{ (£,a,-  e2)  2+  e3} 

e 


da. 
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where  the  change  of  variables 


has  been  used  in  Equation  3-55.  We  can  simplify  sligh 
the  integral  in  question  becomes 


It  is  obvious  that  the  second  integral  yields  a  s 
result  with  a  =  0  and  Ri2  in  place  of  Rix,  namely 


3  8 


Thus  combining  all  these  results,  we  have 


and  since  receiver  decisions  can  be  based  on  ln(l  ) 


Thus  the  most  general  receiver  structure 
generation  of  In  1 ;  is  shown  in  Figure  3-2. 


■UIL> 


UU  ll.TU 


rj  n  »'J  »  .  ■/  .*  .  V  JV’J’STWW'.;  V.  -trrnyr}  ir w-  iru 


FIGURE  3-2.  M-ary  Receiver  Structure  with  Rician 

Fading 


Clearly,  M  structures  of  (generic)  form  as  shown  on 
Figure  3-2  must  be  implemented  and  their  outputs  compared  in 
size  so  that  a  decision  can  be  made  based  on  a  modified 
version  of  the  decision  rule  of  Equation  3-45. 

It  is  obvious  that  if  all  signals  are  equally  likely  and 
Efi  is  independent  of  the  index  i,  the  hardware  inside  the 

dotted  box  shown  in  Figure  3-2  is  unnecessary  so  that  the 
output  from  the  adder  following  the  squaring  devices  is 
sufficient  in  order  to  make  decisions. 

The  receiver  structure  can  be  modified  according  to  the 
constraint  of  a  given  modulation  scheme,  such  as  MFSK  or 
MPSK.  Before  focusing  on  the  performance  of  this  receiver, 
it  is  desirable  to  investigate  what  modifications  must  be 
made  in  the  receiver  structure  of  Figure  3-2  when  5  is  modeled 


l 


* 

i 


.> 

.N 


40 


as  a  random  variable. 


In  order  to  accomplish  this,  we 


backtrack  somewhat  now  and  define 


so  that 


In  (1.) 


where,  as  shown  in  Figure  3-2, 


If  5  is  a  random  variable,  then  1^  as  given  by 
3-63  is  actually  1^18  ,  a  conditional  random  variable 


1, 


f  (8) 
A 


d8 


m  cos (5) 

e 


2 n  I  (m) 


d5 


(3-61) 


(3-62) 

(3-63) 


(3-64) 

Equation 
so  that 


(3-65) 


If  we  examine  separately  the  terms  inside  the  second  set 
of  brackets  in  the  exponential  of  Equation  3-70,  after  some 


manipulations  we  have 


2  2 

r.  +  r , 

1C  IS 


N„ 


on  r.  cos8  -  r.  sin8  2 

ic  +  (X 


2o 


40 


(3-71) 


Thus, 


g. 


1.  =  G  exp  ^  — 


{  -  ♦  frl  } 

L  ^  No  4a  J  J 


1€ 


•jexp| 

-7t 


2a(r,ccos8  -  r.ssin8)  I  0 


m 


cos8 


2<jV 


2rt  I0(m) 


d8 


(3-72) 


Reducing  the  integral  only  yields 

w 

I0(m) 

where 


(3-73) 


(3-74) 


(3-75) 
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g.  +  — ■ 

1  h 


2  2 

V.  +  r< 
ic  Is 


*  fr  * 

40 


m  V<*i> 


In  Ifl(m) 


(3-76) 


The  term  In  Io(m)  cannot  possibly  affect  any  decisions, 


and  the  term  In  Io(qj.)  can  be  approximated  accurately  for 


small  and  large  arguments  of  the  function.  The  generation  of 


In  lj_  ignoring  the  In  I q  ( m )  term  can  be  accomplished  by  the 


system  shown  in  Figure  3-3 


s/2fi(t)C0S[Ci)et+^(t)] 


as 


ijia 


IHn  li 


s/2  fi(  t  )ain[  uet+$<  t  )| 


y/Tp^lnloO 


FIGURE  3-3.  Receiver  Structure  with  Rician  Fading 
and  Random  Specular  Component 


A  significant  simplification  results  if  <<  1.  If  this 


is  the  case,  then 


W  -  T  +  1 


(3-77) 


=>  m  I0(q.)  =  — 


(3-78) 


and 


In  1  ,«g,4 — 

1  1  H 


2  2 

r.  +r.  2 

_ ii  +  _SL_ 

N2  4C* 


1 

4  H 

4 


2a  ri=  T  if 2a 

2a2  No  J  ^2(J2  N0 


-  In  I0(m)  (3-79) 


and  after  some  manipulation. 


In  l.=g.+ 


1  +  SL. 


Mo  4o4N' 


am 


2  2°\ 
r  + - 

ic  « 


+  — CL  2 

^.N2  4^N2 


-  r.  +  r. 

2  ic  is 


am 

—  + - In  I.  (m) 

2  A  0 

4a 


(3-80) 


D .  RECEIVER  PERFORMANCE 

In  order  to  evaluate  the  error  probability  of  the 
receiver,  we  assume  that  sm(t)  is  transmitted  resulting  in  lm' 
taking  on  the  value  1^.  Since 


A 


y  =  —  +  g, 

1  Hi 


(3-81) 


we  see  that 

Pr(c|sm(t) , lm'=Lm}  =  PrUi'ci^,  i*n|  sm(t)  ,  lra'  =Lm) 

In  evaluating  this  probability,  issues  of  independence 
come  into  question.  With  the  assumed  transmitted  signal 
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A 


s  (t) 


+  y  ,  (t )  dt  +  ~^rr 

J  14  2  a 


+ 


+  n  ( t ) 


yi2(t>  dt 


Since 


s  (t) 


r  mk 


(t> 


(3-82) 


(3-83) 


we  have,  after  some  reduction. 


A.  = 

l 


—  [°Ai  *  a  V>]  *  ~  ”, 


_CL_ 

2a 


(3-84) 


where 


•  f 

n.  .  =|  n(t)  y,  . (t)  dt  j=l,2 


(3-85) 


We  see  that  Aj_  is  made  up  of  Gaussian  random  variables 
that  are  squared  and  summed.  We  can  compute  the  p.d.f.  of 

the  component  Gaussian  elements  in  order  to  obtain  the  p.d.f. 
of  Ap  We  let 


b. 

i  ; 


+  -0- 
20 


(3-86) 
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2i 


=  -i-[a,P  •  +  a,£-  1  +  —  n 

71  •  «  2  n  n 


i2 


(3-87) 


Since  a!  and  a2  are  both  Gaussian  random  variables  with  means 
a  and  0  respectively  and  both  with  variance  a2,  we  have 

(3-88) 


e (b. .  (  =  —  at  +  -a- 

11  N.  “  'iJ 
o  2a 


1 


E{b, . )  =  —  aP 

2 1  ^  'mi 


(3-89) 


Also 


b, . -E  {b, . }  =  *““[  (a;-a)  t „  +  a2Pinl  +  -2-  n  , 

•  *  *  *  N_»  ■*  N  “ 

0  0 

v  e<v  ■rI'*,'0,p-‘ +  a^»l  *  r 


(3-90) 


(3-91) 


so  that 


E{  [bu-E{b1.}  ]  [b2.-E{b?i}  ]  }  = 


-  — E  (((»,-«)«,.  ♦  +  n,lH<a,-«»P.i  ♦  a£.  *  n,  J  I 


N„ 


=  —  (02t  P  .  +  C2p.  t  +  n. ,}} 

2  ~im  rmi  ~  im  »:m  li  1 2 

N„ 


(3-92) 


o 


where  the  second  equality  in  Equation  3-92  is  made  possible 
by  the  fact  that  random  variables  a!  and  a2  are  independent . 

Since 

P  .  =  -p.  (3-93) 

rmi  r  im 
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E{nu  ni2 }  =  Eljn(t)  yu(t)  dt  Jn(t)  y  i2  (X)  dx 


‘f 


yn(t)  y,2(t)  dt  *  0 


(3-94) 


it  is  clear  that  b^  and  b2i  are  uncorrelated,  and  since  they 
are  Gaussian  random  variables,  they  are  statistically 
independent.  Furthermore, 


1  2  .2  1  2  q2  1  N0  2 

var{b  . }  =  —  0  t.  +  —  0  B.  +  —  — -  E  =  0 

li  2  nm  2  2  ?  fi  bi 

N  N. 

0  0  0 

.  1  i  J  1  2  r2  1  N.  2 

var(b2.}  =  —  O  +  —  O  5im  +  “  T  ^  '  abi 

No  No  No 


(3-95) 


(3-96) 


Observe  that  var  (b^  )=var  {b2i )  so  that  the  same  symbol  has 
been  used  for  both  variates.  Thus  from  Whalen  [5], 


f  I s  )  = 

A. Is  1  m 


i  V  -J\  “Vi 
— exp '  \  ; —  ■  v 

"bi  l  20-  J  l  J 


(3-97) 


where 


mDi  =  E2{b,.}  +  E2{b21}  =  hr  + 


Since 


if-  It  •>- 


(3-98) 


1. '  =  —  +  g. 

1  H.  1 


(3-99) 
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nmnnmn 


f ,  , ,  (L  ’ Is) 
1.  s  1  m 

i  m 


A  I  s 

l  m 


Li’-  *i 


_1_ 

h 


(3-100) 


1 

MV-  gi)+ 

mi  exp 

2  a. 

bi 

2ob 

bi 

Io 

. 

L.'  >  g. 

i 


(3-101) 


This  obtains  only  a  marginal  p.d.f.  for  li'.  We  want  to 

explore  independence  issues  by  considering  under  similar 
conditions  of  sm(t)  transmitted,  what  correlations  exist  with 
other  variates  bij  and  b2jf  where 


Clearly, 


b,  ,  =  —  (a,^.  +  a-,P  )  +  —  n  + 

15  No  2  ^  No  :1  202 

1  -  1 
b,.  =  —  (aft  .+  a,c.  }  +  —  n.. 


E(b  }  =  —  a$,  +  -SL- 

15  N  >  ,  ? 

o  2(J 


0 

bi  f  E{b!  :}  =  F  {  (V  00  V  a^J»J  +  F  ni’ 


b  -  E{b,.}  =  — {(a,-  a)B  .+  al.  }  +  —  n  . 

2j  2)  N  1  N  ■t 


(3-102) 

(3-103) 

(3-104) 

(3-105) 

(3-106) 

(3-107) 
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Now 


e {  [blfE(bu}]  [bx  -E{bx  }  1 }  = 


-  TEll,a.'  0,1  nnlt<a.-  °>V  *A.+  nl."  <3'1081 

N. 


=  —  [CT2^.  4.  +  P-  +  E{n  n.,}] 

N„ 


(3-109) 


Since 


E { n. ,  n.j}  =  E  Jn (t )  yu<t)  dt  J n (T)  yn  (T)  dx 
>  T» 

No  f  No 

=  TJyii(t)  yn(t)  dt  =  7^ 


(3-110) 


(3-111) 


Thus, 


E  {  [bu-  E{bu}]  [b,.-  E{b1.)l  } 


■il 


N 


o  (4  4.  +  p.  p.  »+  — t 

jm  r  :mr 'm  y  ”  i  ■ 


(3-112) 


which  in  general  does  not  appear  to  be  zero.  Unless  certain 

orthogonality  assumptions  are  made,  it  appears  that  this 

expectation  cannot  be  zero.  Now, 

E  {  [bu-  E  { b , .  }  ]  [b? E(b2;}]  = 


*  —  WIU,-  «>$,.♦  «A.*  "nil  la,-  «»„♦  "„)> 

N. 


(3-113) 


=  —  [a?4.  P  .+  cr2p.  4-  +  E{n.„  n.  )] 

2  rm  j  rim  )2 

N. 


(3-114) 
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Since 


E  { n. ,  nJ2}  =  E  I  Jn  (t )  ytl  (t )  dt  Jn(t)  y  j2  (X)  dt 


(3-115) 


N  f  N 

=  ~  Jyi:(t)  yi2(t)  dt  =  Y  ^ 


(3-116) 


we  have 


1  .  2 


E(  [b.;-  E(b.; )]  [br-  Etb.jl  =  —  [0(tTPT  +  P__4,T)  +yp::l 

N. 


(3-117) 


Similarly,  interchanging  the  indices  i  and  j  in  Equation 
3-117,  we  obtain 


E(  [b.  -  E  {b2 ; }  ]  [bp-  E{b.,m«  —  (<M^Pt>  P^J*  -f  P 

Ph 


(3-113) 


E{  [b?, -  E (b2i }  ]  [b2j-  E{b?.)  ]  i  = 


—  EU;a,-a>pai+  a^im+  n.2][(ar  «>P  +  a£  +  n,?H 


(3-119) 


—  [0?P  P  +  d2^.  %.  +  E { n  n  }  ] 

2  riMrj]  'ui'jm  ;2  ,2 

N„ 


(3-120) 


Since, 


N  r  N. 

h,  n,,>  =  —  y„(t)  dt  =  —  £, 


(3-121) 


we  have 


E{[b?;-  E{b2.}]  [b„-  E{b„}l}»  —CO  <M^+  Wj  +  T- 1 


2i 


2i 


(3-122) 


It  appears  that  unless  some  special  cases  are  considered, 
there  is  some  correlation  amongst  all  the  random  variables 
making  up  1^  '  and  lj'  conditioned  on  the  transmission  of 

s^lt)  • 

Since  the  lj'  cannot  be  demonstrated  to  be  independent,  we 


can  see  that 

Pr  { 1 ;  '  <L~,  l/<Lm, 


.T:“l 


'  <L 


'  <L 


1m’<L  |s  (t ) , 1  ' =L  } 

rl  r*.  T  7* 


can  only  be  approximated  by 
M 

nPr{l.  '  <  L  I  s  (t)  ,  1  '=  L  }  (3-123) 

i  mm  .t  rr. 

i  - 1 , 


where  Equation  3-123  is  an  exact  result  if  and  only  if  the 
random  variables  1 ^ '  can  be  proved  to  be  statistically 


independent .  Therefore 


Pr { c  I  s  (t >  ,  1  '=  L  )  « 

m  fTi  ul 


M 

n  j 

i  -"1, 


1 

g1)+mbi 

-  g  )m: 

-  exp 

K  1 

2 

i . 

a 

u ( L.  -  g  )  dL 


(3-124) 


From  this  we  obtain, 


Pr { c I s  (t) }=]Pr{c|s(t),  1_ ' =  L  }  f 

nl  I  W 


(L  Is  )  dL 


(3-125) 


And  finally 


Pr  { c}  =  Vp  Pr{c|s  (t)  }  (3-126) 

n  m 

While  the  result  is  somewhat  complicated,  simplifications 
can  be  introduced  on  the  terms  that  make  up  the  giver, 
mathematical  expressions. 

E.  SIGNAL  POWERS 

Equation  3-2  gives  an  expression  for  the  signal.  The 
specular  component  of  this  signal  is 

s.  ,  ( t )  a  f .  ( t )  cosi(i)^t+<5  (t )  +  5]  (3-127) 

The  energy  of  this  specular  component  is  found  by 
integrating  Equation  3-127  with  respect  to  time  over  the 
symbol  duration,  Ts,  yielding 

E.  sa2  E,.  (3-123) 

The  diffuse  component  of  the  signal  is  given  by 

s.  D(t)  T  v  f,(t)  cos  [coct+  <(>  (t)  +  8)  (3-129) 

Fixing  v  and  9,  the  energy  of  the  diffuse  component  is  v-  E.;. 

To  find  the  average  energy  of  the  diffuse  component,  wo 
average  this  expression  using  the  p.d.f.  of  the  r.v.  v,  *;  r 
obta in 

E„  =  2  a  E,;  (3  133) 


5  3 


where  for  simplicity  we  will  henceforth  use  the  term  E  to 
represent  signal  energy  in  place  of  since  the  signals 

being  considered  in  the  sequel  have  energies  Efj.  that  are 

independent  of  the  index  "i." 

Thus,  taking  the  ratio  of  specular  energy  (or  power)  to 
diffuse  energy  (or  power) ,  we  obtain 


snrd  2a2 


(3-133) 


1 .  MFSK 

For  M-ary  Frequency  Shift  Keyed  modulation,  the 


signals  take  on  the  mathematical  form 


f, (t)  = 


7! 


i=l , 2, . . . ,M  ;  Ts  =  Tf-  Tfc 


(3-134) 

(3-135) 
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As  a  result  of  this. 


J 


s.  T(t)  s.  T(t)  dt  = 
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7? 


cos 
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“7  +  (l-  7) 


Acot 


co  t  + 
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<•-?) 


Acot 


=  ~  J  cosj  (i-j)  AcotJ  dt  + 


fj 


cos 


2(0  t  +  A(0t 


/  M  \ 

/  .  M  \ 

11-—)  +  Acot 

ID-—) 

\  2  / 

\  2  / 

dt 


The  second  integral  in  Equation  3-136  is 
negligible  regardless  of  value  for  A(o,  so  that 

Tr 

Js.  T  (t)  s  .  T(t)  dt 

sin(i-j)  AcoT  / 2  cos(i-j)  A(0(T  +TJ 

=  E  -  •  - 

(i-j)  A(0Ts/2  2 

We  will  choose 


A(0T 

s 


I 

I 


dt 


(3-136) 


clearly 


(3-137) 


(3-139) 
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where  n  is  an  integer  because  this  choice  guarantees  that  the 
transmitted  signals  are  orthogonal.  Under  these  conditions 
we  have 


lt 

j*  (t)  dt  =  E  =  |  y^2(t)  dt  ;  V  i 


(3-139) 


For  this  modulation  scheme  with  the  assumption  of 
Equation  3-138,  we  have 


*  4 

r  r 

/  .  M  \ 

/  M  \ 

Jr cos 

rfl  S 

Vi)  4at  - 

\m-  t)  Aa>t 

Acotl  dt  =  0  ;  i*m  (3-14  0) 


Tf 

r  f 

/  M  \ 

/  .  M  \ 

i-sin 

Tk  S 

^m-  Acot  - 

I1’  l)  4<at 

(3-141) 


and  after  some  reduction, 


P;  =  E 


sin(m-i)  AcoT  /2  sin(m-i)  Aa)(T  +  T  ) 

_ S  _ D _ I _ 

(m-i)  A(0Ts/2  2 


(3-142) 


Clearly  for  i?tm,  we  will  have  Pim=0f  thus  eliminating 
many  of  the  terms  involving  £j_m  and  in  the  result  for 

receiver  performance  previously  derived. 

Now,  from  the  above  results,  since  Pj_m  =  =  0  for 

i*m.  Equation  3-98  yields 


m,  ,  = 


w 


V  i*m 


(3-143) 


and  Equation  3-95  yields 


2  i. 

0  =  - 

bl  2N„ 


V  i*m 


(3-144) 


In  order  to  obtain  the  probability  of  correct 
(Pr{c}),  we  must  first  evaluate  Pr{c|sm(t),  lm'=  Lm}  as 

provided  in  Equation  3-124,  as  follows. 


L  -g. 

M  m  1 

n  j 


H,l 


2°ol 


^iXmbi 


bi 


u ( x )  dx 
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n 

i*l , i*tn 
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y  e 


2  y2 

~(y  -l) 


I0{y  Ci>  dy 


A (L.-9l 


>aL 


(3-145) 


where  the  change  of  variables 


(3-146) 


has  been  made  in  Equation  3-145. 

The  term  contained  within  the  integral  of  Equation  3- 
145  is  now  clearly  a  Rician  function,  where,  for  convenience, 
we  have  defined 

fiO 

(3-147) 


The  integral  of  a  Rician  in  this  form  is  known  as  a  Marcum  Q- 
Function  and  in  its  general  form  is  defined  as 


oo 


The  probability  of  correct  decisions  for  MFSK  is  now 
obtained  by  applying  Equations  3-125  and  3-126.  Equation  3- 
125  contains  the  p.d.f.  for  the  random  variable  lm'  given 

sm(t)  transmitted,  as  shown  in  Equation  3-101.  The 
parameters  contained  in  this  equation  all  have  the  subscript 
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m,  which  makes  them  different  than  the  parameters  we  used 


earlier  when  the  subscript  was  i  with  i*m.  Reducing 
parameters  in  this  form,  we  have 


=  —  <SNVD 

20 


5T  (SNR»+1) 

“ 0 


SNR, 


2  2 

a  snrd 


<SNRd+1) 


g  =  In  P  -  In (1+SNR  ) 

m  D' 


SNRS 

SNR^ 


where,  again,  the  subscript  "m"  can  be  removed 
equiprobable  signals,  and  Pm  is  simply  M-1 . 

Manipulating  the  above  expressions  into  forms 
useful  for  substitution  into  Equations  3-125  and  3-126, 


SNR, 


bn 


mbn 


bn 


2SNR  (1+SNR  ) 

S  J 


SNR» 


c 


i*i  0 


bn 


=  i  + 


hi  'm 


SNRD 


Thus,  upon  substitution,  the  p.d.f.  of  Equation 


yields 


these 

1-154) 

1-155) 

1-156) 

(-157) 

for 

more 

3-159) 

3-159) 

3-160) 

3-125 
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(3-161) 


Expressing  Equation  3-126  in  its  final  form  in  order 
to  obtain  the  total  probability  of  correct  decisions,  we  have 

2  SNR 


Pr(c)  = 


i 


1  ~  Q 


l  SNRn 


1+SNR  z 

o 


z  G 


z  +2SNR$(l  +  SNiy 


snrd 


2SNRs  (1  +  SNRd) 


SNR^ 


} 


dz  (3-162) 


where  the  substitution 


z  = 


'L*  -  *  T 

®bm 


(3-163) 


has  been  used  in  order  to  obtain  this  result  .  The 
probability  of  error  (P{e})  is  now  obtained  by  subtracting 
from  unity  the  P(c}  shown  in  Equation  3-162. 


2  •  MPSK 


For 

take  on  the 

f, 


M-ary  Phase 
mathematical 


Shift  Keyed 
form 

i  =  l, 2, .  .  ,,M 


modulation , 


T  =  T  -  T 

s  :  c 


the  signals 


(3-16-3) 
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)  (t)  =  -  +  <t> 


1=1,2,  .  .  .,M 


(3-165) 


where  <J>0  is  a  fixed  yet  arbitrary  phase. 

The  signals  T(t)  now  will  clearly  not  be 


orthogonal.  However, 


-  K 


M  (t)  dt  =  y.,(t)  dt  =  E 


and  for  this  modulation  scheme, 


(3-166) 


=  E  cos 


2?C  (i-m) 


(3-167) 


Also, 


(3.  =  E  s 

rim 


f  271  (i-m)  1 

lrl^n 


(3-168) 


Clearly  in  this  case,  none  of  the  terms  involving  and  (3im 
will  cancel  out  as  they  did  in  the  MFSK  case.  After  some 
reductions. 


2n (i-m) 


"bi  “  cTsnrJ  SNRd  +  SNRo  cos|  — 


2  E  SNR3 

CT  =  - 

bi  2N„ 


Thus, 

%  =  2  SNRs 

2  2 

C.  SNR 

bi  D 


(3-169) 


(3-170) 


271  (i-m) 

SNRD  +  cos  ~  + 


(3-171) 


SNR. 


A.  =  - 

im  SNR  +1 


SNR.  +  COS 


[27t  ( i-m)  1  1 

M  J  +  SNR 


(3-176) 


B 


^  J 2(SNRp-H)  (Lm-g) 
SNR. 


C.  = 


2SNR. 


SNR. 


f  f  2rt  ( i-m)  1  1  'l 

\  D  L  M  J  snr2  J 


(3-177) 


(3-178) 


Now,  using  this  notation,  the  p.d.f.  of  Equation 
3—101  reduces  to 


2 

B  -C 

m  mm 


SNR  +1 

f .  , ,  (l  is  )  - - e 

1  '  S  mm  2 

m  m  SNR 


{7?} 


(3-179) 


Expressing  Equation  3-126  in  its  final  form  in  order 
to  obtain  the  total  probability  of  correct  decisions  for  the 
MPSK  modulation  scheme,  we  have 


C.  -A. 

im  im 


Pr  ( c  }  = 


M  M 

Z’.J  n  e 

w*l  i-1 ,  i#m 

—  oo  ' 


{  1  -  Q(A,»  ,  B,o>  } 


B  -  C 

m  mm 


SNRD+l 


SNR. 


n — 1 

/  B  C 

I  m  mm 

V  2 


dB 


(3-180) 


The  final  results  on  receiver  performance  for  MFSK 


and  MPSK  can  now  be  directly  applied  to  a  Rician  fading 


problem.  The  determination  of  the  fading  parameters  that 
identify  the  variables  SNRs  and  SNRD  must  be  determined  either 

through  computer  simulation  or  experimentation  with  a  model 
of  the  reflecting  structures. 

This  last  aspect  of  the  problem  has  been  lert  for 
further  study  as  it  applies  to  the  design  of  the  Space 
Station  digital  communications  links,  since  accurate 
determination  of  the  fading  parameters  was  not  feasible  given 
the  scope  of  this  study. 
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IV.  ANALYSIS  OF  ISI  EFFECTS  ON  M-ARY  RECEIVERS 
WITH  FIXED  MULTIPATH 

In  the  previous  chapter,  it  was  assumed  that  multipath 
interference  resulted  in  signal  fading  effects  causing 
degradation  in  the  performance  of  digital  receivers.  A 
Rician  fading  model  was  used  as  a  generalization  of  Rayleigh 
fading  to  encompass  more  realistic  multipath  scenarios. 

In  this  chapter  analysis  will  be  carried  out  in  order  to 
determine  how  a  single  bit,  delayed  by  multipath  propagation 
and  thus  causing  Intersymbol  Interference  (ISI),  affects 
receiver  performance.  This  different  approach,  no  less 
powerful,  may  be  more  relevant  to  certain  multipath 
scenarios . 

We  analyze  here  the  performance  of  a  receiver  which  is 

optimum  for  for  deciding  which  one  of  M  possible  signals  was 

transmitted,  using  the  following  received  signal  model: 

r  (t )  =  S.(t)  +  n(t)  t0St<Ts+t0  ;  i=l,2,...,M  (4-1) 

where  Ts  is  the  symbol  duration,  s^(t)  is  a  deterministic 
signal,  and  n(t)  is  once  again  the  AWGN  of  PSD  level  Nq/2 
Watts/Hz . 

The  receiver  which  is  optimum  for  this  problem  is  shown 
in  F igure  4 - 1 . 
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FIGURE  4-1.  Optimum  M-ary  Receiver  Configuration,  No 

ISI 


It  is  now  intended  to  derive  the  performance  of  this 
receiver  when  multipath  is  present  so  as  to  introduce 
transmission  delays  not  to  exceed  one  symbol  time.  The 
direct-ray  path  signal  will  be  assumed  to  have  zero  delay  so 
all  multipath  contributions  will  have  delays  referenced  tc 

the  delay  of  the  direct-ray  path  signal.  For  mathematical 
simplicity,  let  t0=0  and  assume  that  Sj(t)  is  transmitted. 


Then,  in  the  absence  of  multipath, 
r (t)  =  s( (t)  +  n{t) 


0  <  t  <  T. 


Consider  the  effect  of  multipath  now.  We  introduce  the 
following  notation 


S  (t  ;(X,N,T)  = 


S  .  ( t  -  T  ) 


[♦It] 


j  =  1,2,  ...  ,M  a„>0  x_>0 

where  there  are  N  secondary  ray  paths  available  at  the 
present  symbol  time  during  which  the  signal  travels  from  the 
transmitter  to  the  receiver,  and  a„  and  are  the  strengths 

and  delays  respectively  of  the  signal  which  arrives  via  these 
secondary  paths. 

If  Sj(t)  is  transmitted,  with  multipath  present, 

r(t)  =  s.(t)  +  s  .  (t ; (X,  N,X  )  +  s^ft+T;  j},L,rp  r  n(t)  ( 4  -  ) 

where  0<t<T,  j  =  l,2,...,M,  and  k=  1 , 2 ,  .  .  .  ,  M .  Equation  4-'-, 
assumes  that  in  the  previous  symbol  time,  the  symbol  sk(t) 

was  transmitted  during  which  there  were  L  secondary  ray 
paths,  with  strengths  PL,  P2,  •••,  Pl  and  corresponding  delays 
T|  1 ,  tj2 ,  .  .  .  ,  til.  Note  that  L  does  not  necessarily  equal  N ,  ar 
does  not  necessarily  equal  pn,  and  Xn  does  not  necessarily 
equal  Tln.  Note  that  since  for  the  previously  transmitted 
symbol  to  arrive  at  the  receiver  at  the  same  time  as  the 
currently  transmitted  symbol,  the  latter  symbol  would  have  to 
have  taken  a  different  path,  so  the  attenuating  and  dela\ 
factors  corresponding  to  the  current  and  previous  symbol 
all  likelihood  may  not  be  equal. 

In  all  analysis,  since  X  represents  the  delay  r 
current  bit  and  T|  represents  the  delay  of  the  previous  tit, 
we  will  assume 

Xi<x2<...<xn  ;  0<xr.<Tn 
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Tii-Ti2-’  •  ^11 


T<Tl1<2Tn 


To  consider  a  very  specific  case,  let  us  assume  that  we 
are  dealing  with  a  modulation  scheme  in  which  all  signals 
have  equal  energy  and  equal  probability  of  being  transmitted, 
namely 

T 

Ei  =  |  s*(t)  dt  =  E  Vi  (4-5) 

o 

P,  =  -  Vi  (4-6) 

1  M 

where  T  will  henceforth  be  used  to  express  the  symbol 
duration,  Ts,  and  the  signal  set  is  orthogonal,  that  is 

T 

Js.(t)  SJ(t)  dt  =  E  Si5  V  i,  j  (4-7) 

o 

Thus,  assume  that  Sj(t)  is  transmitted,  sk(t)  was 
previously  transmitted,  and  define  (see  Figure  4-1) 

T 

li  sjr(t)  s .  ( t )  dt  Vi  (4-8) 

o 

Under  the  assumed  conditions,  the  probability  of  a 

correct  decision,  denoted  Pr(c  I  j , k }  becomes 

Pr{c  I  j , k }  =  Pr{l  >1.;  V  i*j  |  j , k }  (4-9) 

and 

r(t)  =  s  (t)  +  s  .  (t ;  a,N,  X)  +  sv  (t+T;p,  L,rj)  +  n(t)  (4-10) 

0  <  t  <  T 


where 


Thus 


+  J  n  (t)  s;  (t)  dt 

o 


J 


1  N 

s.(t)  s,(t)  dt  +  I  2^  an  s,(t-T_)  s.(t)  dt 
o'  1  0  n=1 


n  n*l 


s.  (t+T-  tl  )  s.  (t)  dt  +  n. 

n  k  *n  l  i 


(4-11) 


Hi5} 


n(t)  s.(t)  dt 


(4-12) 


1.  =  E  8  t  + 

1  1  j 


z«j 

r.-l  „ 


S(t-X)  S(t)  dt 


ih\ 


+  /  13  I  (t+T-  Tln)  s .  (t)  dt  +  n 

w 

n*i 


(4-13) 


1,  = 


T 

jl 


s,  (t)  +s.  <t;a,N,T)  +sk  (t+T 


;P,L,n)] 


S  (t)  dt 


In  general,  the  remaining  integrals  in  this  expression 
will  not  be  zero,  even  if  i*j,  and  i*k .  For  example,  for 
XF3K  modulation, 


cos  (2Jt  f  t) 


S  (t) 


0<t<T 


( 4  -  .  4 1 


and 


J 


s. (t)  s,(t)  dt 


■tJ 


COS(27Cf.t)  COS(27tf,t)  dt 


sin  [2rt  (f .  -  f  )  T]  sin[2ji(f.+  f  )  T  ] 

=  E  - - - — - - -  +  £  3 


2JC(f.~  f .)  T 


2Jl(f.+  f.)T 


(4-15) 


If  and  ( f  ^  +  f  j )  T  each  equal  some  integer,  then 

indeed 


1  ** 


(t)  s.(t)  dt  =  E  8.. 


(4-16) 


The  frequencies  fj_  and  fj  can  be  chosen  to  satisfy  Equation  4- 

16  and  this  will  be  assumed  to  be  the  case  in  the  sequel. 
Furthermore, 


I 


■if 


s  (t-t  )s.(t)dt  =  | —  cos  [2JCf  .  (t-t  )  )  cos  [2jif .  t  ]  dt  (4-17) 

T  J  n  • 


which  after  some  manipulation,  reduces  to 


{  sin  (27tf  T  Tn  *  ]  -  sin[2Jtf<Ttn' ]  } 


27CT  (  f ;  +f  .) 

{  sin  [  2Jtf  ,T  Tn '  ]  -  sin(2nf.TXn']  }  s  E^tt.1)  (4-18) 


2jcT  (  f  -f  ) 


s 


il  so, 


s  (t+T-TlJ  s  (t)  dt 


cos  [2icf<  (t  +  T-  T|  )  ]  cos[2rtf  t]  dt 
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27CT(fk  +  f . ) 


2  7tT  ( f  k  -  f , ) 


{  sin[27CT(fk+f ,  tln'  )  ]  -  sin[2jtfkT(Tln’-l)  ]  } 
.  {  sin[27CT(fk-f.  Tln’)  ]  +  sin[27CfkT(T|n,-l)  ]  } 


s  E  *) 


(4-19) 


where 


0<T  ’<1 

n 


(4-20) 


1<T1  '<2 


(4-21) 


are  normalized  delays. 


Thus, 


u 

■i  -  E  6i,+X“„  VV  tlk<V>  +  n, 


(4-22) 


Clearly,  lt  is  a  conditional  Gaussian  r.v.  for  which  its 

mean  and  variance  can  be  computed.  Under  the  assumption  of 
Sj(t)  presently  and  sk(t)  previously  transmitted.  Equation 
4-22  specifies  1^  I  j,k,  and  using  an  overbar  to  denote 
statistical  expectation,  we  have 


1,  I  j,  k  =  E 


8n  ♦£p„Ci„<V> 

n  - 1  n  -  i 


(4-23) 


1  i  I  j,  k  -  1;  I  j,k  1J  j,k  -  1.1  j,k 


n  n 


=1/ 


N  Nq 

— -  8(t-  x)  s  (t)  s  (t)  dt  =  —  e  5 

o  l  m  o  im 


(4-24) 


From  this  we  observe  that 


No 

var ( 1.  I  j,  k }  =  —  E 


(4-25) 


Thus,  the  random  variables  1^1 j,k  are  statistically 
independent,  so 

M 

Pr  { c  |  j ,  k }  =  J~Jpr{l  .>1.  |  j ,  k  } 


i=l, i*  j 


“  M 

■  I  n  fi  ij.k'v 


dL. 


(4-26) 


where  Lj  is  the  value  assumed  to  be  taken  on  by  lj.  Since 

L, 


Pr  { 1.<L  .  |  j 


if  k)  =  f  -r  „ 1  ,  r  ■ , ■;  exp 

-oo  V27tN0E/2 


(x-  1. I j, k) 


2  N  E/2 

o 


dx 


=  erf 


n- 1 


(4-27) 


where 


erf # (v) 


■J 


■J~ZK 


u 

~  T 

S  du 


(4-28) 


Thus , 
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Pr {c | j, k}  = 


M 

J  n. 

-co  l-i.w 


L  -  E  i  (n,  i,  j,  k) 


rN0E/2 


r27tNQE/2 


2 

3  {V  E  r<n'i'^k)  }  dL 

N0E  ^  3  *  : 


(4-29) 


where 


Pi 

r  (n,  i,  j ,  k)  =5,,  +  X0tn^ij(xn,)  +  XM,<V>  (4_30) 


We  obtain  finally, 


P  r  { c  }  = 


tt 

j  1  k«l 


Pr  { c  I  j  ,  k  }  Pr  { s  .  (t )  ;  s<(t)} 


(4-31) 


where  Pr{Sj(t)  ;  sk(t)}  is  the  probability  that  Sj(t)  is 
presently  transmitted  and  s^Ct)  was  previously  transmitted. 

Due  to  the  independence  of  transmissions,  the  probability 
that  sk(t)  was  transmitted  during  the  previous  symbol 

interval  is  the  same  as  the  probability  that  it  will  be 
transmitted  during  the  current  symbol  interval.  Since  the 
symbols  have  been  assumed  to  have  equal  probability, 


Pr  {  s  .  (t )  )  =  Pr  {  s  (t)  }  =  — 
3  *  M 


we  have 


The  probability  of  error  (P{e})  is  now  obtained  by 
subtracting  the  probability  of  correct  decision  shown  in 
Equation  4-30  from  unity. 

In  order  to  check  this  result,  assume  now  that  no 
multipath  is  present,  so  that  an  =  Pn=0  V  n.  Then  the  above 

expression  correctly  yields  the  well-known  probability  of 
correct  reception  of  M-ary  orthogonal  signals,  namely 


Again,  our  derived  result  is  only  valid  for  the  specific 
cases  that  satisfy  Equations  4-5  through  4-7.  This  means 
that  nonorthogonal  modulation  schemes  such  as  MPSK,  QAM, 
etc.,  are  not  covered  by  these  results.  To  include  the  more 
generalized  cases  of  nonorthogonal  modulation  schemes  would 
require  a  similar  yet  mathematically  more  complex  approach 
and  therefore  has  not  been  attempted  here. 
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A.  DESCRIPTION  OF  THE  RELEVANT  PORTIONS  OF  THE  SPACE 

STATION  COMMUNICATIONS  SYSTEM 

Figures  5-1  and  5-2  show  diagrams  of  the  designed 
physical  configuration  of  the  National  Aeronautics  and  Space 
Administration  (NASA)  Space  Station  (SS)  as  of  3  May,  1986. 
This  spacecraft  will  represent  a  milestone  in  the  evolution 
of  space  research  and  manufacturing. 

Incorporated  into  the  Space  Station  design  is  an 
elaborate  multiple  access  communication  system.  This  system 
will  provide  for  constant  communications  with  a  number  of 
earth  stations  by  maintaining  links  through  NASA's  Tracking 
and  Data  Relay  Satellite  System  (TDRSS)  .  The  system  design 
incorporates  the  ability  to  simultaneously  communicate  with  a 
number  of  space  platforms,  including  Co-orbiters/Free-Flyers 
(FF),  the  Orbital  Transfer  Vehicle  (OTV) ,  the  Orbital 
Maneuvering  Vehicle  (OMV) ,  the  Space  Shuttle  Orbiters,  and 
astronauts  conducting  Extra-Vehicular  Activity  (EVA) . 

Astronauts  conducting  EVA  will  be  suited  in  a  newly- 
designed  spacesuit,  designated  the  Extravehicular  Maneuvering 
Unit  (EMU) .  The  EMU  will  have  communications  capabilities  by 
maintaining  links  for  voice  and  data,  as  well  as  video  links 
for  television  cameras  and  for  relaying  a  Heads-Up-Display 
(HUD)  projection  from  the  Space  Station  into  the  EMU  helmet. 


FIGURE  5-1 


NASA  Space  Station  Baseline  Configuration 


The  Space  Station  is  currently  in  the  process  of  full- 
scale  development.  Correspondence  with  Mr.  Sid  Novosad  of 
NASA's  Johnson  Space  Center  [6]  is  the  sole  reference  for  the 
technical  descriptions  contained  in  this  chapter. 

As  will  later  become  apparent,  the  only  communications 
link  of  interest  insofar  as  the  development  of  this  chapter 
is  concerned  is  the  Space  Stat ion-to-EVA  astronaut  forward 
link.  There  are  two  system  configurations  currently  under 
consideration  for  this  link.  The  pertinent  characteristics 
of  these  two  configurations  are  outlined  in  Table  5-1. 

TABLE  5-1.  SS  TO  EMU  FORWARD  LINK 


Modulation 


BPSK 


(Alternate) 

Frequency  13.7  GHz 

Bit  Rate,  Data  40  KBPS 

Bit  Rate,  Video  400  KBPS 

Link  Margin  3  dB 


DCPSK 

QFSK 

14.7  GHz 
100  KBPS 
22  MBPS 
6.1  dB 


Consider  first  the  communication  links  required  for  the 
Space  Station  to  transmit  to  the  Co-orbiters/Free-Flyers , 
Orbital  Transfer  Vehicle,  Orbital  Maneuvering  Vehicle,  and 
Shuttle  Orbiters.  Since  the  transmitting  and  receiving 
platforms  are  in  free-space,  devoid  of  reflectors,  the  only 
potential  source  for  multipath  is  the  structure  of  the 
receiving  platform  itself.  In  order  to  determine  if  the 
receiving  platforms  are  "large"  enough  to  provide  sufficient 
reflected  signal  travel  delay  so  as  to  cause  multipath 
propagation  loss,  we  must  first  determine  how  far  the  RF 
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energy  will  travel  during  one  bit  interval. 


This  is 


obviously  a  function  of  bit  rate.  Assuming  propagation  at 
the  speed  of  light  in  free-space,  we  make  use  of  the  bit 
rates  under  consideration  in  order  to  generate  Table  5-2. 


TABLE  5-2.  DISTANCE  TRAVELLED  IN  ONE  BIT  PERIOD 


Bit  Rate.  Rb  (BPS )  Bit  Period.  T  Distance, 


22 

M 

45.45  ns 

13  . 

400 

K 

2.5  (is 

750 

100 

K 

10  (is 

3000 

40 

K 

2  5  (is 

7500 

Ht— Lml 


The  proposed  receiving  platforms  for  Space  Station  are 
substantially  smaller  than  the  corresponding  signal 
propagation  distances  for  all  except  the  highest  bit  rate 
case.  In  addition,  all  platforms  are  generally  devoid  of 
components  in  their  structure  that  would  cause  reflection  of 
a  substantial  amount  of  RF  energy  directed  towards  the 
receive  antenna.  Thus,  this  situation  prevents  multipath 
from  being  a  serious  problem,  even  for  the  22  MEPS  case. 

Nevertheless,  the  Space  Station  is  a  relatively  large 
structure,  in  relation  to  the  13  meters  transit  width 
required  to  produce  ISI  given  a  22  MBPS  signaling  rate.  The 
dual  keels,  the  living  and  working  modules  (denoted  the 
Habitability  Modules),  the  solar  panels,  and  the  solar 
dynamic  reflectors  are  all  reflectors  of  RF  energy.  This 
means  that  multipath  is  a  concern  when  the  receiving  antenna 


is  physically  close  to  or  within  the  confines  of  the  Space 
Station  structure. 

Some  comments  must  be  made  here  regarding  the  use  of 
"space  diversity"  as  a  multipath  countermeasure.  In  more 
traditional  atmospheric  fading  situations,  the  propagation 
ray  paths  experience  refractive  "bending"  within  the 
atmosphere.  If  all  the  RF  energy  from  a  signal  is  bent  away 
from  a  particular  area  in  the  atmosphere,  that  area  is  called 
a  "shadow  zone,"  because  the  signal  is  not  easily  received 
within  that  area.  This  effect  can  be  overcome  by  placing 
receiving  antennae  at  different  locations  so  that  if  one 
antenna  falls  within  the  shadow  zone,  others  may  still  be 
able  to  receive  the  signal.  This  technique  is  called 
"spatial  diversity." 

When  a  radio  signal  is  transmitted  to  the  Space  Staticr., 
the  conductive  elements  which  make  up  its  frame  can  reflect, 
diffract,  absorb,  and  re-radiate  the  RF  energy  so  that  the 
ray  paths  of  the  signal  are  bent.  Thus,  with  the  Space 
Station,  shadow  zones  may  exist,  but  not  because  of  the  i:.  : 
of  refractive  ray  trace  bending  that  occurs  :  r.  * 
atmosphere.  Nevertheless,  spatial  diversity  t  echn  i  :  ;*■  s  e  ■■ 
still  useful  in  overcoming  these  signal,  losses.  S  h'A.'A  . 
planning  to  use  a  similar  technique  on  the  Spare  .‘a*  :  • 
Current  plans  are  to  distribute  receiving  ar.tenr.  i- 
extreme  of  its  structural  periphery  s-~  :r.*er  •  .  :•  : 

before  they  are  bent  by  the  st  ruct  urai  e  ler  . 


SO 


This  spatial  diversity  technique  also  allows  us  to  rule 
out  severe  multipath  loss  in  most  cases  when  the  Space 
Station  is  the  receiving  platform.  If  an  antenna  intercepts 
the  incoming  signal  before  it  is  reflected  by  the  Space 
Station  structure,  then  the  antenna  cannot  experience  the 
multipath  which  is  caused  by  that  structure.  An  onboard 
computer  will  process  the  received  signals  in  real-time  in 
order  to  determine  which  antenna  is  receiving  the  signal  with 
the  least  amount  of  interference,  and  will  link  that  antenna 
to  the  receivers  directly.  Thus,  multipath  is  not  expected 
to  be  a  serious  problem  when  the  receiving  platform  is 
located  on  the  Space  Station. 

This  leaves  one  case  yet  to  be  considered.  If  an 
astronaut  involved  in  Extravehicular  Activity  is  positioned 
within  the  periphery  of  the  Space  Station  structure,  the 
radio  transmissions  have  ample  opportunity  to  suffer  multiple 
reflections  before  they  arrive  at  the  astronaut's  own  receive 
iroer.na.  The  direct  path  will  undoubtedly  exist,  yet  there 
an  re  a  plethora  of  alternate  paths  reflected  off  all  the 
:  ."tat  ion  extremities  previously  listed.  Furthermore,  m 

a  small  structure,  basically  a  "space  suit”,  ref 
:  :  v.  i :  r.  g  'he  opportunity  of  he  i  r.g  able  to  make  use  of  spate 
r  y  as  a  countermeasure  for  multipath.  In  addition, 

•  t  ‘  :.»*  a-'f»  station  sh^rt  -  ra.oae  antenna  and  the 
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accurately  determine  whether  multipath  will  in  fact  seriously 
interfere  with  the  astronaut's  ability  to  communicate  with 
the  mother  ship. 

A  complete  analysis  cannot  be  performed  here.  But  it  may 
be  possible  to  judge  qualitatively  whether  this  subject 
deserves  further  attention  by  evaluating  previously  obtained 
results  on  receiver  probability  of  error  when  ISI  effects  due 
to  multipath  propagation  for  conditions  typical  of  what  can 
be  expected  for  Space  Station  designs.  That  is  the  intent  of 
this  chapter. 

B.  ANALYSIS  OF  THE  POTENTIAL  FOR  SPACE  STATION  ISI 
1 .  Scenario 

In  order  to  evaluate  the  probability  of  error 
expression  given  by  Equations  4-31  and  4-32,  it  is  necessary 
to  set  up  an  arbitrary  geometry  for  which  the  requisite 
assumptions  that  lead  to  multipath-induced  ISI  are  met.  Due 
to  analytical  constraints,  however,  the  geometry  must  be 
arranged  so  that  the  only  bits  arriving  at  the  receiver  at 
time  t  are  the  current  bit  and  the  previously  transmitted  bit 
only.  A  complete  analysis  would  require  a  detailed  survey  to 
determine  exactly  how  many  paths  satisfy  this  criterion,  and 
to  establish  the  expected  signal  strengths  and  delays 
associated  with  each  path.  However,  if  it  can  be  shown  that 
a  single  reflected  path  can  interfere  sufficiently  with  the 
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quality  of  communications  as  to  effectively  nullify 


receiver,  then  a  more  detailed  analysis  may  be  unnecessary. 

Assuming  that  the  very  large  solar  panels  of  t 
Space  Station  are  the  primary  RF  reflectors,  and  assurr.i 
that  a  22MBPS  data  link  is  in  use  between  the  Space  Star: 
omni  antenna  and  an  EVA  astronaut.  Figure  5-3  descric^s 
geometry  where  the  "kth"  bit  arrives  via  the  direct  pa' 
while  at  the  same  time,  the  previous  "k-lth"  bit  arrived 


the  reflected  path. 
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time,  the  previous  bit  "k-1"  arrives  via  the  reflected  path, 
and  after  a  slight  delay,  the  "kth"  bit  follows. 


In  this  scenario  it  is  assumed  that  the  solar  panel 
:s  facing  away  from  the  astronaut,  so  that  its  aluminum 
rackplate  faces  him/her.  Since  aluminum  is  a  good  conductor, 
it  is  possible  that  approximately  100%  of  the  incident  14." 
i  H .’  carrier  wave  is  reflected  by  the  solar  panel, 
r  -  r  “  hern'  re,  since  *•  he  space  .oss  difference  between  the 
lire:-,  and  reflected  paths  is  negligible,  the  reflected  path 

■;jr,al  has  approximately  the  same  electromagnetic  strength  as 
*  r.e  iiiect  path  signal.  Thus,  the  an  as  well  as  the  [i. 
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Shift  Keying) ,  QFSK  (Quaternary  Frequency  Shift  Keying) ,  and 


DCPSK  (Differentially  Coherent  Phase  Shift  Keying) .  The 
analysis,  however,  does  not  apply  to  a  differential 
modulation  scheme,  and  although  it  may  be  possible  to 
extrapolate  results  for  differential  schemes  from  the 
analysis,  it  would  be  quite  difficult  to  do  so  in  the  general 
case,  therefore  such  schemes  are  not  considered. 

Although  BPSK  is  being  considered  as  one  of  the 
potential  modulation  schemes  for  Space  Station  digital 
communication  links,  that  scheme  is  not  an  orthogonal 
signaling  technique,  so  that  the  general  results  obtained  in 
this  chapter  do  not  apply  to  this  specific  scheme.  However, 
since  the  performance  of  most  systems  employing  BPSK 
modulation  enjoy  a  3  dB  SNR  improvement  over  the  performance 
of  similar  systems  utilizing  BFSK  modulation,  which  is  an 
orthogonal  signalin  scheme,  the  results  involving  the  latter 
scheme  can  be  used  as  an  indicator  of  a  similar  system 
employing  BPSK  modulation  provided  that  approximately  3  dB 
SNR  adjustments  are  made. 

The  receiver  probability  of  error  has  been  evaluated 
min?  *  he  general  results  of  Chapter  4,  for  the  cases  -f 
•  r  irc.m :  ns  i'.n  via  BFSK  and  QFSK  modulated  signals,  as  \ 
■  "  :  *  -n  of  SNR  .  In  this  chapter,  SNR  is  defined  as 
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where  E  is  the  signal  energy  and  N0  is  the  PSD  level  of 
Additive  White  Gaussian  Noise  (AWGN)  in  Watts/Hz.  The  signal 
transmission  and  reflection  geometry  used  in  order  to  obtain 
these  results  is  based  on  the  assumptions  described  in  the 
previous  section.  The  signal  strengths  traveling  along  the 
reflected  paths,  namely  a  and  p,  were  assumed  to  be  unity, 

indicating  the  fact  that  we  are  assuming  that  the  strength  of 
the  signals  traveling  along  the  reflected  path  equals  that  of 
the  signals  traveling  along  the  direct  path.  The  delays,  t 
and  ti  ,  were  calculated  to  be  0.348  given  the  assumed 
propagation  geometry. 

Numerical  performance  evaluations  were  then  performed 
in  which  the  signal  traveling  along  the  reflected  path  had  a 
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FIGURE  5-5  OFSK  With  Delevs 


Keeping  now  the  signal  propagation  delay  constant  j' 


0.1,  Figures  5-6  and  5-7  display  the  corresponding  results 


for  signal  transmission  via  BFSK  and  QFSK  when  a  and  (i  vary 


over  the  range  from  zero  to  one.  These  plots  show  -ha*  t  : 


this  relatively  small  signal  propagation  delay,  systems 


BFSK  modulation  display  a  marginal  amour.*  <.:  f  *-• ;  f  rm.an 


superiority  over  systems  using  QFSK  modulation. 
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FIGUR*  5-7.  QFSK  With  D«l«y«0 . 1 


In  general,  for  this  limited  scenario,  ISI  analysis  shows 
that  : 

•At  a  22MBPS  data  rate,  degradation  can  result  even  from 
just  one  multipath  source  even  given  the  fact  that  in  this 
analysis  only  one  bit  was  interfering  at  any  given  time. 

•ISI  is  worse  when  the  delay  is  long  (see  Figures  5-6 
through  5-9) . 

•ISI  gets  worse  as  the  reflected  signal  strength  gets 
stronger . 

•In  the  worst  cases  of  ISI,  the  differences  between  QFSK 
and  BFSK  are  too  small  to  be  of  concern. 

Although  just  one  reflected  path  was  analyzed,  the 
probability  of  error  was  shown  to  have  quickly  degraded  to 
unacceptable  levels.  In  light  of  the  fact  that  Space  Station 
will  have  much  more  than  one  multipath  reflector  and  will  not 
:;e  limited  to  just  one  bit  interference,  a  nominal  design 
specification  link  margin  of  6  dB  may  not  reduce  the 
probability  of  error  to  acceptable  levels.  Therefore,  this 
is  a  problem  that  deserves  more  careful  scrutiny  before  the 
actual  hardware  choices  are  made. 
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Chapter  2  provided  a  brief  overview  of  the  root  causes  of 
multipath  propagation.  The  general  theory  pertaining  to  how 
multipath  propagation  degrades  digital  receiver  performance 
was  described.  In  addition,  a  description  of  how  the 
physical  structure  relates  to  the  analytical  parameters  was 
provided . 

Chapter  3  then  dealt  with  multipath  propagation  in  the 
form  of  a  digital  communication  system  experiencing  Rician 
fading.  Analytical  expressions  of  receiver  probability  of 
error  for  the  specific  cases  of  MFSK  and  MPSK  modulation 
scnemes  were  obtained. 

Chapter  4  treated  the  multipath  propagation  phenomenon  in 
the  form  of  single-bit  Intersymbol  Interference  (  I S I )  .  A 
exact  result  for  receiver  probability  of  error  was  obtain*' 
for  M-ary  orthogonal  signaling  schemes  where  the  M  symhc 
have  equal  probabilities. 

In  Chapter  5,  multipath  propaaat  ion  analysis  was  aprli** 
to  NASA's  Space  Station.  It  was  demonstrated  that  in  i 
least  one  communication  link  scenario,  In tersymh 
I.nterference  may  cause  unacceptable  receiver  performance  in  a 
link  between  the  Space  Station  and  an  Extravehicular  Activity 

It  must  he  emphasized  that  Chapter  a  t 


( EVA)  ast  ronaut . 


only  a  cursory  but  conservative  look  at  the  actual  conditions 
to  be  present  when  Space  Station  is  finally  built.  Actual 
receiver  performance  could  be  somewhat  worse  than  the 
performance  obtained  in  this  study. 

A  further  complication  is  that  in  the  analysis  of 
Chapter  5,  the  signal  delays  and  signal  strengths  were 
treated  as  deterministic  quantities.  In  actuality,  since 
there  will  be  no  way  of  predetermining  the  physical  locaticr. 
of  the  astronaut's  receiver  at  any  instant  in  time,  it  would 
be  more  accurate  to  treat  those  parameters  as  random, 

processes . 

Since  the  phase  of  an  incoming  signal  is  determined  by  its 

propagation  delay,  and  since  the  propagation  delay  in  the 

case  of  the  Space  Stat  ion-to-ast ronaut  link  is  most  likely  * 
be  observed  as  a  random  process,  the  delay  itself  should  : 
treated  as  a  random  process.  In  any  situation  where 
incoming  signal  phase  is  random,  it  is  generally  unwise  * 
encode  the  message  as  a  function  of  that  phase.  This  n  v. 

preclude  the  use  of  Phase  Shift  Keyed  signals  a.  : : 

appropriate  modulation  scheme  for  Space  Station.  In 
s  i  *•  uat  ion  such  as  the  one  described  in  Chapter  c,  a  V :  y 

Sh  i  t  -  Keyed  modulation  scheme  is  usually  a  ref  ter  ;h  i 
since  the  frequency  of  the  incoming  signal  will  he  re  la*.  : 
unaffected  by  the  signal's  propagation  over  such  a  oh  :  • 
d i  s t  an.ee  . 
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